i ADA049176 


PROJECT  SQUID 

TECHNICAL  REPORT  UC-1-PU 


MECHANISMS  OF  OECAY 
OF  LAMINAR  ANO 
TURBULENT  VORTICES 


MAHINDER  S.  UBEROI 
UNIVERSITY  OF  COLORADO 
BOULDER,  COLORADO  80309 


PROJECT  SQUID  HEADQUARTERS 
CHAFFEE  HALL 
PURDUE  UNIVERSITY 
WEST  LAFAYETTE,  INDIANA  47907 


DECEMBER  1977 


D D c 

JAN  25  1978 


’irn: 


Project  SQUID  it  a oooparative  program  of  basic  rasaarch  relating  to  Jat  Propul- 
sion. it  is  sponsored  by  the  Office  of  Naval  Research  and  it  administered  by 
Purdue  University  through  Contract  N00014-7S-C-1143,  NR-098-038. 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


/i-efpu]j)-  / 


Technical  llepart^^UC-l-PU 


PROJECT  SQUID* 


A COOPERATIVE  PROGRAM  OF  FUNDAMENTAL  RESEARCH 
AS  RELATED  TO  JET  PROPULSION. 

OFFICE  OF  NAVAL  RESEARCH*  DEPARTMENT  OF  THE  NAVY 
CONTRACT/Tra0014-75-C-n43i  NR-098-038 


MECHANISMS  OF  D 
AND  TURBULE 


ECAY  OF  LAMINAR 
NT  VORTICES^ 


M Mahinder  S./Uberoi  / 
University  of  coiorad^o 
Boulder,  Colorado  80309 


PROJECT  SQUID  HEADQUARTERS 
CHAFFEE  HALL 
(PURDUE 

WE?r  LAFAYETTE,  INDIANA 


This  document  has  been  approved  for  public  release  and  sale; 
its  distribution  is  unlimited. 


3 b // 


Mechanisms  of  decay  of  laminar  and  turbulent  vortices 


By  MAHINDER  S.  UBEROI 

University  of  Colorado,  Boulder,  Colorado  80309 

Dynamics  of  an  infinitely  long  one-dimensional  vortex  and  a swirl  are  com- 
pared with  dynamics  of  a semi-inf inltely  long  trailing  vortex  and  trailing 
swirl.  With  increasing  distance,  the  change  in  axial  velocity  difference  between 
the  core  of  the  tralling-vortex  and  the  surrounding  region  causes  radial 
convection  and  some  associated  axial  convection  of  angular  momentum.  In  laminar 
or  turbulent  trailing  vortices,  this  is  the  dominant  mechanism  for  decrease  in 
velocities  of  swirl  in  the  core  and  corresponding  growth  of  the  core.  Axial 
velocity  difference  between  the  core  of  the  tralling-vortex  and  the  surrounding 
region  is  necessary  for  the  sustenance  of  turbulence  in  the  vortex  core.  A 
theory  of  turbulent  tralling-vortex  is  developed  on  the  basis  of  these  mechanisms 
and  the  results  are  compared  with  experimental  observations. 


1.  Introduction 


There  are  a variety  of  theories  and  views  about  unconfined  and  seml- 
infinitely  or  infinitely  long  turbulent  swirls  and  vortices.  These  theories 
are  neither  confirmed  nor  refuted  by  experimental  investigations,  since  very 
few  experiments  exist  where  the  effects  of  initial  conditions  and  extraneous 
influences  have  been  minimized.  We  critically  examine  the  assumptions  under- 
lying these  theories  and  show  that  many  of  these  are  far-fetched.  We  find 
correct  physical  processes.  The  results  of  analysis  based  on  these  physical 
processes  are  compared  with  the  most  recently  available  experiments. 

There  are  four  distinct  types  of  flows  which  are  relevant  here.  An  Infinitely 
long  LINE  SWIRL  with  a single  velocity  component  Ug(r,t)  , where  r Is  tlie 
radial  distance  from  the  swirl  axis  z , subscript  0 refers  to  the  angular 
coordinate,  and  t Is  the  time.  Its  angular  momentum  pM  per  unit  axial 
distance  is  constant  and  Independent  of  time,  where 

M = 2tt^/  Ug(r,t)r^dr  , (1) 

and  p is  the  density  of  the  fluid.  A prime,  as  u'  , will  be  used  to  denote 
the  fluctuating  part  of  a quantity.  Otherwise,  the  symbol  represents  either  a 
quantity  in  a laminar  flow  or  its  mean  value  in  a turbulent  flow. 

The  second  flow  is  an  infinitely  long  LINE  VORTEX  with  one  velocity 
component  Ug(r,t)  such  that  Ug  = r^/2iir  for  large  r with  zero  velocity 
at  r “ 0.  Its  angular  momentum  per  unit  distance  along  the  swirl  axis  is  in- 
finite and  its  rate  of  change  of  angular  momentum  is  finite  depending  on  vls- 

2 

coslty.  It  is  further  assumed  that  3r  /3r>0  where  T « 2Trru„ . 

0 


1 
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The  third  flow  Is  a TRAILING  SWIRL  with  three  velocity  components  u^(r,z), 

u (r,z),  and  u (r,z),  where  z is  the  distance  along  the  axis  of  the  trailing 
o r 

swirl  measured  from  its  origin.  It  may  be  produced  by  rotating  vanes  with  the 
axis  of  rotation  parallel  to  the  prevailing  uniform  steady  flow  u^.  The  vanes 
may  add  to  or  absorb  some  axial  momentum  of  the  prevailing  flow.  Thus  we  have 
a swirl  which  may  have  a coaxial  jet  or  a wake.  As  r u^  = u^  = 0, 

and  u^  = u^.  The  flux  of  angular  momentum 

A = 2Trp^/“u^(r,z)  Ug(r,z)  r^  dr  (2) 

is  Independent  of  z. 

The  fourth  flow  is  a TRAILING  VORTEX  produced  at  the  tip  of  a semi-infinite 

lifting  wing  in  the  presence  of  a prevailing  mean  flow.  The  trailing  vortex 

has  some  similarity  to  a trailing  swirl  having  three  velocity  components,  with 

the  important  difference  that  as  r->-«>,u  = 0,  u.  = r /2iir  , u = u 

where  T is  the  total  circulation.  The  flux  of  angular  momentum  is  infinite, 
0 

but  its  change 


dA/dz  = 2np  J "^(u  u.)  r dr 


3z'“z  e'  ‘ 

2 

is  finite  and  depends  on  z . It  is  further  assumed  that  3r  /3r  > 0. 

Initially,  in  all  these  flows  vorticity  is  assumed  to  be  confined  to  a 
region  of  small  radius  and  moments  of  all  orders  of  any  vorticity  component 
with  respect  to  z-axis  exist.  The  cores  of  these  flows  are  regions  where 
most  of  vorticity  is  located.  This  paper  is  concerned  with  the  growth  of  these 
cores  or  the  spread  of  initially  concentrated  vorticity  due  to  viscous  and 


turbulent  processes. 


w 


r ^ 


In  analyzing  the  third  and  fourth  flows,  it  is  invariably  assumed  that  u^ 

may  be  neglected,  u = u , and  z is  replaced  by  tu  , thus  reducing  the 
z o o 

trailing  swirl  to  the  line  swirl  and  the  trailing  vortex  to  a line  vortex.  We 
indicate  below  and  subsequently  show  in  detail  that  this  approximation  is 
invalid. 

f 

As  a trailing  swirl  or  a trailing  vortex  develops,  swirl  velocities 
Ug  decrease  with  the  downstream  distance  z.  Since  the  pressure  at  large  r 
is  constant,  near  the  axis  of  a swirl  or  a vortex  this  leads  to  3p/8z  > 0 
where  p is  the  pressure.  This  implies  divergence  of  the  cores  of  these 
flows  and  an  axial  velocity  difference  between  the  cores  of  these  flows  and 
the  surrounding  regions  that  have  three  important  effects  on  the  dynamics  of 
these  trailing  flows  which  are  absent  in  a line  swirl  and  a line  vortex. 

(1)  Mnear  and  non-linear  stability  analyses  (Uberol  et  al.  1972;  Narain  and 
• Uberoi  1973)  show  tliat  a difference  of  axial  velocities  between  the  core  and 

I the  surroundings  destabilizes  swirling  Clows  which  otherwise  could  be  stable, 

j (2)  There  may  be  significant  and  sometimes  dominant  radial  and  associated  axial 

I 

I convection  of  angular  momentum.  (3)  The  range  of  downstream  distances  over 

p 

j which  dynamic  self-similarity  exists  may  be  limited.  The  importance  of  these 

[ effects  decreases  with  decreasing  rate  of  spread  of  these  flows.  However,  in 

i 

j the  study  of  trailing  flows  with  swirl  the  emphasis  is  on  their  rates  of  growth 
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laminar  flow,  a part  of  which  may  become  turbulent.  In  swirling  flows  the 
basic  laminar  flow  may  enhance,  diminish,  or  even  quench  turbulence  in  its 
interior.  These,  sometimes  strong,  stabilizing  or  destabilizing  effects  must 
be  considered  when  postulating  turbulent  stresses  in  these  flows. 

The  equation  of  motion  for  a swirl  or  a line  vortex  is 


3 2 3 3 3 ./, 


(4) 


where  v is  the  kinematic  viscosity. 

A swirl  of  finite  M diffuses  out  due  to  viscosity  and  shares  its 
angular  momentum  with  the  surrounding  fluid  which  is  set  into  motion.  Its 


Reynolds  number  is  (M/t)  /v.  A known  self-similar  swirl  is 


''  “ = 2^  ^ ^ (-r2/4vt). 


(5) 


In  a line  vortex  u„  = F / 2wr  for  r ■>•  «>  and  the  rate  of  change 

u O 


of  angular  momentum. 


It 


r u„dr  “ -2vl’ 

0 o 


(6) 


This  rate  is  Independent  of  the  detailed  distribution  of  the  vorticity 


or  u in  the  interior  of  the  vortex.  If  we  assume  finite  v then  there  are 

y 


stresses  but  no  net  force  on  a fluid  element  in  the  potential  flow  surrounding 
the  core  where  most  of  the  vorticity  resides.  The  angular  momentum  is  lost 
from  the  Interior  of  the  vortex  through  the  potential  flow  region  to  the  region 
r ->  ® . The  core  over  which  most,  say  95%,  of  the  total  vorticity  or  F is 


distributed  grows  and  so  does  its  angular  momentum.  This  Is  due  to  infinite 


■ V*-5 
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angular  momentum  surrounding  any  finite  though  growing  interior  region.  The 
Reynolds  number  is  • 

A known  self-similar  solution  for  the  line  vortex  is 

r = 2TrrUg  = ~ exp(-r^/4vt)  ^ (7) 

We  may  combine  (5)  and  (7)  to  get  a vortex-swirl  combination 

^ = ^0  ^ 2^  - 1 } exp  (-r^/4vt)31  (8) 

0 

In  a line  vortex  (7)  shows  that  the  circulation  increases  monotonlcally 
with  radius  reaching  a constant  value  . In  a swirl  total  vorticity  is  zero 
and  r increases  and  then  decreases  to  zero  for  large  r.  In  a line-vortex-swirl 
example  given  by  (8)  F overshoots  F^  , then  decreases  to  F^  for  r ► . 

This  overshoot  decreases  and  becomes  relatively  insignificant  with  time  and  may 
be  considered  as  a decaying  'Initial'  disturbance. 

This  is  a special  case  of  a general  result.  If  the  Initial  vorticity  is 
axisymmetric  and  is  distributed  over  a finite  area  around  the  origin  then 
its  subsequent  distribution  can  be  found  using  (4).  After  a long  time  vorticity 
and  velocity  distributions  will  appear  as  if  the  total  vorticity  were  originally 
at  the  origin,  any  initial  vorticity  distribution  of  opposite  signs  and  zero 
total  value  having  no  significant  influence.  If,  however,  the  total  vorticity 
is  zero,  then  the  initial  distribution  of  the  vorticity  determines  the  sub- 
sequent state  of  the  vorticity  and  the  velocity. 

We  may  look  at  the  situation  from  the  point  of  view  of  dynamics.  The 
angular  momentum  associated  with  a finite  total  vorticity  is  infinite  and  it 
is  finite  for  a distribution  of  vorticity  of  opposite  signs  such  that  the  total 
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vorticity  or  is  zero.  As  time  progresses  the  former  will  dominate  the  latter. 
We  are  assuming  initially  concentrated  vorticity  near  the  axis. 

These  well  known  results  are  presented  to  contrast  some  properties  of 
laminar  flows  with  corresponding  properties  of  the  flows  with  the  same  overall 
parameters  and  in  turbulent  state.  For  example,  in  a laminar  line-vortex- 
swirl  the  dynamics  of  any  region  of  overshoot  where  F > F^  become  in  time 
unimportant  to  the  dynamics  of  the  main  vortex.  While  in  a turbulent  line- 
vortex-swiri  Covindraju  and  Saffman  (1971)  assert  that  at  high  Reynolds  numbers 
the  overshoot,  F > F^ , is  the  main  growth  mechanism  of  a turbulent  vortex. 

The  virtual  origin  of  the  time  t for  a line-swirl  and  line-vortex  may 
not  be  the  same.  This  becomes  insignificant  as  t becomes  large.  However, 
for  both  small  and  large  t we  define  a line-vortex-swirl  as  a flow  which 

behaves  like  a vortex  for  sufficiently  large  r and  for  some  finite  range  of 

2 

r,  dr  /dr  < 0.  This  slight  generalization  is  necessary  for  the  purpose  at 
hand,  which  is  to  study  the  stability  and  turbulence  in  such  flows. 

An  Important  criterion  based  on  analysis  assuming  invlscid  flow  for 
stability  of  swirling  flows  with  one  velocity  component  u^^  (r,t)  is,  (Rayleigh 
1916),  (Chandrasekhar  1961), 

dF^  /dr  > 0.  (9) 

A line  vortex  is  stable  at  all  times.  A swirl  is  unstable.  A llne- 
vortex-swirl  is  unstable  and  its  interior  may  become  turbulent  at  a sufficiently 
high  Reynolds  number.  However,  as  the  time  progresses  the  flow  will  tend  to 
stabilize  and  production  of  turbulence  will  diminish  towards  zero  and  the 
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turbulence  already  created  will  decay  due  tu  viscosity.  In  practice  flows 
approximating  a line-vortex  are  used  to  stabilize  unstable  fluid  conf Igura- 
tions  such  as  vortex  stabilized  electric  arc,  (Chow  and  Iberoi  1972).  Experi- 
ments lend  support  to  the  above  criterion  (9)  without  regard  to  any  limitations 
of  the  cited  stability  theories. 

We  have  conducted  a simple  experiment  to  check  the  stability  considera- 
tions. Water  was  injected  tangentially  all  along  the  inner  wall  of  a transparent 
vertical  cylindrical  vessel  of  30  cm  diameter,  50  cm  long,  with  a central 
drain  in  its  flat  bottom  and  nearly  full  of  water.  After  the  vortex  was  set 
up  the  drainage  was  reduced  in  steps  after  each  repetition  of  the  experiment 
described  below.  There  was  no  drainage  during  the  last  experiment.  The  fluid 
in  the  center  was  made  turbulent  by  stirring  it  randomly  or  by  spinning  a 0.6 
cm  diameter  rod  spanning  the  entire  length  of  the  cylinder  in  the  direction  of 
the  main  vortex  or  opposing  it.  The  turbulent  flow  was  made  visible  by  painting 
the  stirrer  or  the  rod  with  water  soluble  ink.  Care  was  taken  to  limit  the 
disturbance  to  a short  time  period,  so  as  to  confine  the  initial  disturbance 
to  a cylindrical  region  of  about  3 cm  diameter,  which  may  be  considered  small 

in  size  compared  with  the  size  of  the  main  vortex.  The  Reynolds  number  f /v 

3 

of  the  vortex  was  about  10  . In  every  case  the  initial  turbulence  decayed  and 
was  not  sustained  at  the  expense  of  the  energy  of  the  relatively  slowly  changing 
main  motion  of  the  vortex. 

4 

In  another  experiment  at  F /v  = 7.8  x 10  we  have  made  detailed 

0 

velocity  measurements  in  a trailing  vortex  when  it  approached  a line  vortex 

for  a short  distance,  l.e.  u u . The  existing  turbulence  in  the  vortex 

z o 
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almost  disappeared  (Singh  and  Uberoi  1976).  This  further  substantiates  our 
claim  that  turbulence  cannot  be  sustained  in  a line  vortex. 


3.  Turbulent  Line-Swirl,  Line- Vortex,  and  Their  Combination 


A line-swirl  is  unstable  according  to  the  criterion  (9)  based  on  the 
assumption  of  an  inviscid  fluid.  At  sufficiently  high  Reynolds  numbers  it 
becomes  unstable,  which  is  consistent  with  experience.  We  have  recently  de- 
veloped the  structure  of  a turbulent  swirl  based  on  plausible  assumptions, 
(Uberoi,  1977a). 

A line-vortex  is  stable  unless  it  has  a swirl  superimposed  on  it; 

2 

i.e.,  there  is  a finite  radial  region  where  dT  /dr  < 0.  Turbulence  will 
develop  in  the  nature  of  an  initial  disturbance  and  decay  at  a faster  rate 
than  the  asymptotic  rate  of  growth  of  the  vortex  core.  However,  there  are 
several  theories  of  sustained  turbulence  in  a line-vortex.  They  assume  that 
there  exists  three-dimensional  turbulence.  However,  there  is  only  one  mean 
velocity  component  u^  (r,t).  The  equation  governing  a turbulent  line-swirl, 
line-vortex,  or  their  combination  Is 


Dt 


8r 


2 

r 


"o 


+ vr 


3r 


(10) 


where  - pu  'u„'  is  the  turbulent  stress, 
r 0 

We  note  that  the  equation  (6)  for  the  rate  of  change  of  angular  momentum 
of  a turbulent  line-vortex  or  llne-vortex-swirl  is  the  same  as  for  the  laminar 
case.  Since  any  turbulence  is  confined  to  the  vortex  core,  it  does  not  affect 
the  ultimate  transfer  of  angular  momentum  through  the  potential  flow  surrounding 
the  vortex  core. 
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Squirt-  (1965)  was  the  first  to  conslclt-r  a vortex  with  a turbulent 
core.  In  effect,  he  assumed  that  "turbulent"  kinematic  viscosity, 

V = af  (11) 

t 0 

where  a is  a constant.  The  solution  is  given  by  (7)  with  v replaced  by 

of  . Since  the  turbulence  is  confined  to  a finite  radius,  v should  vanish 

o t 

as  r ->•  “>.  Squire's  assumption  allows  far  too  much  angular  momentum  to 
escape  to  infinity  which  is  determined  by  (6). 

We  may  try  to  save  Squire's  solution  by  stipulating  that  it  is  only 
valid  for  a vortex  in  the  presence  of  uniform  atmospheric  turbulence  where  a 
constant  may  be  used.  However,  there  are  serious  difficulties  with  this 

artifice.  A constant  turbulent  kinematic  viscosity  due  to  atmospheric  turb>s 
lence  has  nothing  to  do  with  which  is  associated  with  the  vortex.  Let 

the  atmospheric  turbulence  be  strong  enough  to  Interact  with  the  vortex.  In 
the  potential  part  of  tlie  vortex,  due  to  spatially  varying  rate  of  strain, 
the  interaction  would  vary  spatially  and  with  time.  Consequently  v^  cannot 
be  assumed  constant.  Further,  we  cannot  assume  that  the  entire  vortex  inter- 
acts significantly  with  the  atmospheric  turbulence,  and  the  outer  flow  is  still 

potential  with  u.  = F /iTfv  for  all  time, 
y 0 

Hoffman  and  Joubert  (1963)  considered  radial  transfer  of  angular  momentum 
in  the  turbulent  core  of  a line-vortex.  Using  certain  assumptions  they  concluded 
that  the  circulation  varies  logarithmically  with  radial  distance  in  the  region 
of  maximum  u^  . They  failed  to  show  how  the  momentum  is  transferred  to  large 
radial  distances;  this  transfer  is  determined  by  viscosity  and  is  given  by  (6). 


w 
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The  total  rate  of  change  of  angular  momentum  is  due  to  viscosity  and 
cannot  exceed  that  given  by  (6).  If  we  insist  on  the  growth  of  the  core, 


i.e.,  decrease  of  its  swirl  velocities  faster  than  that  caused  by  viscosity, 

then  the  outer  flow  must  speed  up,  since  the  total  angular  momentum  must  be 

conserved  except  for  a small  loss  due  to  viscosity.  It  follows  that  the 

circulation  in  the  region  of  potential  flow  where  the  flow  speeds  up  must 

exceed  or  over-shoot  F . Various  elaborate  theories  have  been  developed  to 

o 

"prove"  the  existence  of  a circulation  over-shoot.  In  fact,  the  over-shoot 
is  a direct  consequence  of  the  insistence  mentioned  above,  which  may  take 
many  different  forms. 

2 

Covindaraju  and  Saffman  (1971)  assume  that  r u ' u„ ' and  ru„  are 

r 6 0 

c 

functions  of  r/t*  and  (4)  becomes  a total  differential  equation.  The  insistence 
is  contained  in  this  functional  dependence;  i.e.,  the  maximum  value  of  u^ 
must  decrease  and  the  core  size  must  increase  faster  than  would  be  caused  by 
viscosity. 

Macagno  and  Macagno  (1975)  assume  that  if  e is  the  mean  rate  of  strain 
then  the  turbulent  kinematic  viscosity 

Vt  = a + 4^  ^^2) 

where  a and  6 may  vary  with  space  and  time  but  are  taken  to  be  constant  in 
their  analysis.  The  quantity  a here  is  not  related  to  that  in  (11).  The 
above  is  supposed  to  Include  vortex  generated  and  atmospheric  turbulence.  In 
accordance  with  the  above  discussion  of  Squire's  work  a must  equal  v . 

[ 

i 

I The  equation  (12)  allows  turbulent  stresses  and  hence  production  of  turbulence 

I in  the  outer  potential  flow  where  e is  finite. 

i 
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Theories  discussed  here  and  other  such  theories  were  based  on  the 
unjustifiable  belief  that  a line-vortex  is  equivalent  to  a trailing  vortex 
where  t is  replaced  by  • The  results  of  the  analyses  were  compared 

with  experimental  observation  of  trailing  vortices. 

Measured  data  on  trailing  vortices  have  been  fitted  to  Squire's  solution 
for  a line-vortex,  although  the  value  of  the  constant  a varies  from  case 
to  case  (Rose  and  Dee  1963;  McCormick  et  al  1968).  Owen  (1970)  has  given  an 
explanation  for  the  variations  of  a with  / v . This  does  not  remove 
the  fundamental  objections  raised  above. 

4.  Laminar  Trailing  Swirl,  Trailing  Vortex,  and  Their  Instabilities 


The  equation  governing  u in  these  flows  is 

0 


L- 

3z 


r ^ ^ 2^3.  2, 

'^0’^  + 3?  > 


3 3 /0, 


(13) 


where  u^  is  the  constant  prevailing  velocity  along  the  swirl  axis.  In  the 
core  of  these  flows  the  swirl  velocity  u^  decreases  with  downstream  distance 
7.  and  therefore  > 0.  The  flow  in  the  core  diverges,  causing  significant 

radial  and  some  associated  axial  convection  of  angular  momentum.  We  illustrate 
this  by  examining  the  total  rate  of  change  of  angular  momentum  in  a trailing 
vortex.  Integrating  (13),  we  have 

2iru  / — u„r  dr  = -2ir  — (u  u.)r  dr-  furl  r-+'“>.r  - 2vr 

oo-*  3z  0 3z  z 0 L r J 0 0 

convection  diffusion 

(U) 
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Batchelor  (1964)  has  calculated  axial  flow  in  a trailing  vortex,  neglecting 
radial  flow.  He  assumes  that  all  velocities  are  small  compared  with  u^  , so 
that  (13)  becomes 


u 

o 


3z 


3 33 

IT  > 


(15) 


which  replaces  a trailing  vortex  with  a line-vortex.  The  velocity  u^  is 
given  by  (7)  with  t replaced  by  • Using  this  Ug  pressure  is  cal- 

culated from  the  approximate  equation 


2 , ^ 


3r 


(16) 


Axial  velocity  u is  calculated  using  this  pressure  and  the  equation 
z 


3 3p  1 / ^ 

pu—  u =--:^+pv(- 
o 3z  z 3z 


3r 


2 r 3r  ^ “z 


(17) 


The  result  is 


where 


r zu 

o „ o 

u = - :: — icn  

z ...  2 V 


32it  vz 


r Lu 

exp  (-0  + — ^2  ~ 8vi“  exp(-f.), 

32tt  z 


(18) 


C = u^  r / 4vz 


(19) 


The  function  Q2  and  the  constant  L are  not  of  interest  here.  We  have 
calculated  the  radial  and  associated  axial  convection  as 


r .(u  r) 
o r r ->■  “ 


-2  / u Tr  dr  = r /16n  zu 
O'*  3z  z 0 0 


(20) 
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The  Importance  of  the  neglected  convection  of  angular  momentum  can  be 
expressed  as  the  ratio 

convection/diffusion  = ^^^rdr  + r^*(u^r)^  (21) 

= (r^/v)  (r^/647rzu^) 

= (c/400z) 

where  we  have  assumed  that  the  trailing  vortex  is  generated  by  a semi-infinite 

wing  of  chord  c and  T = cu  /2  , Since  the  vortex  is  essentially  a high 

o o 

Reynolds  number  phenomenon  1 at'd  z/c  must  be  large  to  make  this 

ratio  much  smaller  tlian  unity.  Tn  general  radial  and  associated  axial  con- 
vection of  angular  momentum  cannot  be  neglected. 

Moore  and  Saffman  (1973)  have  calculated  axial  velocity  In  the  (!ore 


of  a vortex  for  wliich  at  z = o,  u^  = br  where  (1  Is  a constant  and 
0 < n < 1.  They  also  neglected  radial  and  associated  axial  convection  of 
angular  momentum.  Using  their  axial  velocity,  we  find  that  the  requirement  for 
these  neglects  is  that 

( ^ )"  « 1 (22) 

0 

They  were  concerned  with  axial  velocity  during  vortex  sheet  roll  up  near  the 
wing  or  small  z where  the  flow  is  essentially  three-dimensional  and  u^  can 
not  be  neglected.  The  condition  (22)  may  be  satisfied  at  large  Z , but  then 
vortex  is  rolled  up  and  corresponds  to  Batchelor's  case. 

The  present  discussion  shows  that  for  those  distances  of  Interest  from 
the  origin  where  the  vortex  is  significantly  changing  the  dominant  mechanism 


13 


for  a decrease  in  the  swirl  velocities  in  the  core  of  a laminar  trailing 

vortex  is  radial  and  associated  axial  convection  of  angular  momentum. 

Another  important  effect  of  axial  flow  or  difference  in  u between 

z 

the  core  and  the  surroundings  is  that  a flow  becomes  unstable  which  was 

otherwise  stable  according  to  criterion  (9),  (Uberol  et  al.  1972). 

A trailing  vortex  is  further  destabilized  when  a trail ing-swirl  is 

added  to  it  such  that  there  is  a finite  radial  distance  r for  which 
2 

dr  /dr  < 0 , in  the  same  manner  as  for  a line-vortex. 


5.  Turbulent  Trailing  Vortex 


On  dimensional  grounds  we  may  write  the  functional  dependence  of  the 

circulation  T ( - 2Trru-)  as 

0 

r/r^  ■ r /r,  , . (23) 


All  "free"  turbulent  flows,  l.e.,  without  constraining  rigid  boundaries,  are 
independent  of  Reynolds  number.  There  has  been  some  speculation  (McCormick 
et  al.  1968)  that  the  structure  of  a turbulent  trailing  vortex  may  depend  on 
its  Reynolds  number  • 

The  equation  governing  u using  slender  vortex  core  approximation 

0 


is 


8 . ^ , 2 3 

^ = 37 


(-U  u 

r 0 


-u  'u. 


+ vr 


3r  r ^ 


(24) 


Integrating  this  equation  we  have 
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(25) 


u /^T^rrdr  = - / 7“  u I'rdr  - (u  r)  F - 2vr 
0 O''  Dz  O'*  3z  z r r >■  “ o o 


«ooH 

/ |-  u (F-r  ) rdi 

o'  dZ  Z 


- 2vF 


(26) 


convection  diffusion 

where  we  have  assumed  that  turbulence  in  the  trailing  vortex  vanishes  as 
r ->■  “ and  we  have  made  use  of  the  continuity  equation  in  deriving  (26).  It 
is  well  known  that  in  all  flows  without  constraining  boundaries  the  turbulent 
fluid  is  separated  from  non-turbulent  fluid  by  a sharp  irregular  boundary. 
Velocity  fluctuations  decay  very  rapidly  as  we  move  from  turbulent  into  non- 
turbulent  fluid,  Further,  measurements,  Singh  (1974),  Uberol  (1974),  In  a 


2 h 2 h 

trailing  vortex  show  that  both  r(u  ' ) and  r(u  ' ) vanish  as  r 

u r 

2 

Hence  r u'  u'  ^ o as  r ->■  “>  . We  are  concerned  here  with  the  spread  of 
r 9 

turbulence  which  is  initially  concentrated  near  the  axis  of  the  trailing  vortex. 

In  all  cases  of  turbulent  flow  thus  far  experimentally  Investigated,  the 

Irregular  front  separating  the  turbulent  from  non-turbulent  fluid  propagates 

2 

at  a finite  rate  rather  than  diffusing  to  infinity.  The  vanishing  of  r 

is  consistent  with  all  known  experimental  facts  about  the  trailing  vortex  and 

2 

the  turbulent  flows  in  general.  Assuming  that  r is  finite  as  r -*■  ® 

would  lead  us  to  the  same  difficulty  as  Squire's  work  (1965).  Angular  momentum 
would  escape  to  infinity  due  to  turbulent  stresses  far  in  excess  of  that  allowed 
by  laminar  viscosity  of  the  fluid. 

It  follows  from  (26)  that  the  radial  and  associated  axial  convection  of 
angular  momentum  are  Important  if  the  turbulent  vortex  grows  faster  than  the 
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W 


I 


laminar  vortex  and  there  Is  no  overshoot  or  P . The  experiments 

determine  the  magnitude  of  the  relative  importance  of  convection  versus 
diffusion. 

The  velocities  u and  u in  a turbulf*nt  trailing  vortex  behind 
6 z 

an  airfoil  have  been  measured  by  Singh  (1974)  and  Uberoi  (1974)  at 

4 

r /v  = 2.1  X 10  . Unfortunately,  the  convection  of  angular  momentum  cannot 

o 

be  accurately  calculated  from  measured  u^  . We  have  calculated  the  first 
and  the  last  terms  in  (25)  using  measured  Ug  and  thus  determined  that 
diffusion  is  about  a percent  or  so  of  tlie  convection  of  the  angular  momentum. 
Therefore,  the  dynamics  of  turbulent  trailing  vortex  are  independent  of 
F^^/v  at  least  at  the  Reynolds  number  of  the  experiment  and  above.  In  the 
literature  the  effect  of  slowly  decaying  and  different  initial  conditions 
in  different  experiments  may  have  been  confused  with  the  effect  of  Reynolds 
number  on  the  structure  of  these  vortices.  The  terms  involving  v may  be 
neglected  in  (23)  through  (26)  (Uberoi  1977b);  thus 


r/r  = Y(ru  /r  , zu  /r  ) 

o oo  oo 


u — u„  r = -;r-  r (-u  u.  - u u ')  - u u r 
0 9z  0 3r  r 6 r 6 9z  z 6 


In  order  to  proceed  further  we  may  write  the  equations  governing  u^ 
and  u^  and  assume  enough  relations  among  the  independent  variables  so  that 
their  number  equals  the  number  of  equations.  Instead  we  propose  an  elemental 
tlieory  which  incorporates  the  mechanism  of  vortex  changes  discussed  above. 

We  look  for  a solution  such  that  is  a function  of  a single  variable 


* m. 


r u 2 ^ 


(29) 


0 o 

We  assume  that  the  total  radial  and  associated  axial  convection  of  angular 
momentum 


o/“  [I?"'  Ve] 


zu 


m 2 

• r 


(30) 


The  terms  on  the  right  hand  side  of  (28)  are  only  significant  in  the  core 
where  F - is  significantly  different  from  zero.  The  sign  of  these 
terms  should  not  depend  on  the  sign  of  F - F^  and  they  should  have  proper 
dependence  on  r as  r ->■  0.  On  these  bases  and  in  light  of  discussion  of 
the  physical  phenomena  we  propose  that 


_3_ 

ar 


(u  u-  + u'u')+  — r u u 
r 6 r 0 3z  z 0 


n a ( 


o 

zu 


( r - r ) 

O 


— exp  (bn) 


(31) 


o 

where  a and  b are  constants.  The  factor  exp(bn)  is  in  recognition 
of  the  fact  that  the  turbulent  core  at  any  z is  of  finite  size  and  the 
expression  (31)  should  rapidly  approach  zero  as  we  go  from  turbulent  core  to 
tlic  outside  non-turbulent  fluid. 

Using  (31)  the  governing  equation  (28)  becomes 


u Fr  = -n  a ( — - )*"  ( F - F ) 

o 3z  zu  0 

o 


2 

exp  (bn) 


(32) 


and  for  m = 1 - n we  have 


±1 

dn 


2 

a(l-Y)  exp(bn) 


(33) 
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The  solution  is 


Y = 1 - 1/(1  + ^ [exp(bn)  - i]  ) (34) 

The  comparison  of  this  expression  with  observations,  Singh  (1974),  Uberoi 
(1974)  is  shown  in  figure  1,  where  a = 150  and  b = 10. 

6.  Final  Stages  in  Vortex  Decay 

If  in  (30)  m > 0 then  relative  to  diffusion  the  importance  of  convection 
of  angular  momentum  decreases  as  z Independent  of  the  theory  proposed 

here,  let  us  assume  that  the  convection  is  negllgable  compared  with  the  dif- 
fusion and  the  core  size  continues  to  increase  at  least  at  the  rate  given  by 
the  diffusion  as  z With  increasing  core  size  and  no  radial  and  associated 

axial  convection  of  angular  momentum,  the  axial  velocity  becomes  unimportant. 

We  clafm  that  under  these  conditions  tlie  flow  becomes  stable  ami  no  sustained 
turbulence  Is  possible.  See  section  2 above. 

In  studies  of  laminar  and  turbulent  line  vortices  self-similarity  is 
often  assumed  (Squire  1965,  Govlndraju  and  Saffman  1971).  Consider  the  fol- 
lowing form 

r/r^  = Y(r^ct")  = y(0  (35) 

Integrating  (10)  and  using  (35)  we  have 

nP  t*^'^ 

It  follows  that  n = 1 and  c ' v.  Using  (35)  the  expression  for  the  rate  of 
change  of  kinetic  energy  of  mean  motion  becomes 
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(37) 
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_9 

9t 
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dr 


2 

-r  /8TTt 

o 


2 

It  may  be  shovm  that  I /Sirt  is  the  rate  of  viscous  dissipation  of  kinetic 

o 

energy  of  a laminar  line  vortex  given  by  (7) . In  a line  vortex  with  a tur- 
bulent core,  the  rate  of  decrease  of  kinetic  energy  should  exceed  that  in  a 
laminar  line  vortex  with  the  same  over-all  parameters.  Hence,  even  if  a tur- 
bulent line  vortex  exists  it  cannot  have  the  form  given  by  (35) . We  may  say 
that  the  virtual  origins  of  t are  different  for  laminar  and  turbulent 
line  vortices.  However,  this  becomes  unimportant  as  t 


7.  Discussion 


We  have  shown  that  all  previous  theories  of  decay  of  trailing  vortices 
are  deficient  or  Just  plain  wrong  since  they  are  based  on  untenable  hypotheses. 
We  find  the  dominant  physical  phenomena  of  radial  and  associated  axial  convec- 
tion of  angular  momentum  and  the  role  of  axial  velocity  in  sustaining  turbu- 
lence in  the  vortex  core. 

A theory  for  a turbulent  trailing  vortex  is  presented  which  satisfies 
the  requirement  of  turbulent  theories.  We  propose  an  expression  for  the 
dominant  physical  phenomena  which  is  consistent  with  the  basic  equations  and 
the  results  agree  with  observations. 

Reasonably  complete  discussion  is  presented  because  of  confusion  in  this 
field  and  to  provide  suggestions  for  further  experimental  work  in  turbulent 
vortices.  It  is  necessary  to  conduct  more  extensive  experiments  to  accurately 
determine  the  values  for  m and  n.  Once  we  have  accurate  measurements  of 
u^  then  we  may  use  the  following  relation  to  independently  determine  m , 
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'I 


\i 


, V-  , 


from  (31)  we  have 


o o 


In  the  past  there  have  been  no  guiding  physical  processes  or  theories 
which  could  help  evaluate  various  devices  and  methods  for  amelioration  of 
vortex  wake  problem  and  its  influence  on  the  operation  of  airplanes  which  may 
interact  with  trailing  vortices  of  other  airplanes.  It  is  hoped  that  present 
theory  and  discussion  of  physical  phenomena  will  provide  such  guidance. 

This  research  was  initiated  during  a University  of  Colorado  Faculty 
Fellowship. 
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We  have  used  Batchelor's  work  (1964)  to  estimate  the  radial  convection 
and  the  associated  axial  convection  of  angular  momentum.  His  work  lias  been 
criticized  by  Tam  (1973),  which  in  turn  has  been  criticized  by  Herron  (1974) 
We  can  estimate  the  terms  neglected  in  the  linear  theory  directly  from  the 
equations. 

The  linearized  equation  of  motion  for  u^  is 


u 
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p 3z 
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A1 


In  the  linear  theory  it  is  further  assumed  that  the  pressure  Is  given  by  the 
equation 
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p Or  “ r 


where 


2r„er  . r^d-e-") 


and 
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n = u r /4vz. 
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This  leads  to 
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p 3z 
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8vz 
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A3 


where 


P(n)  = /" 


A4 


For  the  sources  given  by  A3  the  diffusion  equation  A1  has  no  self-similar  or 
asymptotic  solution  (see  below) . Let  us  assume  that  at  z 

o 


f 


u (r  ,z  ) = 0 

Z 0 


Due  to  the  pressure  gradient  u^  will,  of  course,  change.  This  is  a 

reasonable  initial  condition  since  the  main  interest  is  the  change  u due  to 

z 

the  prescribed  pressure  gradient. 

The  radial  convection  and  associated  axial  convection  of  angular  momentum 


at  z^  is  given  by  (14),  thus 


^ [ T"  u (r-r  ) r dr  = /“ 
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where  we  have  made  use  of  A1  - A5  and  the  value  of  the  definite  integral 
is  approximately  1/2.  The  ratio  (see  eq.  21) 


convection/diffusion  = (F  /v)(r  /bAirz  u ) 

0 0 0 0 


which  agrees  approximately  with  (21). 

In  oruer  to  examine  the  general  solution  of  A1  we  consider  the  change  in 
the  Integrated  velocity  which  is,  of  course,  determined  by  the  sources  given 
liy  Al.  Int>'gr/i(  lug  Al  and  using  A1  and  A4  for  large.  ii,  we  h.ave 


00 


or 
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A9 


where  we  have  assumed  that  there  is  no  velocity  deficit  at  z^.  The  integrated 
velocity  deficit  given  by  A9  continues  to  Increase  with  z and  can  be  made 
to  have  any  value  by  suitable  choice  of  z^.  A particular  self-similar 
solution,  namely 
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8vz 


e J P(t)e*'dt 
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AlO 


was  proposed  by  Tam  (1973)  and  gives  Infinite  Integrated  axial  velocity 
deficit  as  it  is  should  in  view  of  our  discussion. 

For  the  sources  given  by  A3  the  only  properly  posed  problem  is  to  prescribe 
an  initial  u^(r,z^)  and  A1  will  determine  the  subsequent  development  of  u^. 
There  is  no  solution  which  is  either  self-similar,  or  independent  of 
u^(r,z^)and  z^  (depending  only  on  the  integrated  u^(r,z^)). 
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Dynamics  of  an  infinitely  long  one-dimensional  vortex  and  a swirl  are  com- 
pared with  dynamics  of  semi -infinitely  long  trailing  vortex  and  trailing 
swirl.  With  increasing  distance,  the  change  in  axial  velocity  difference  be- 
tween the  core  of  the  trailing-vortex  and  the  surrounding  region  causes  radial 
convection  and  some  associated  axial  convection  of  angular  momentum.  In 
laminar  or  turbulent  trailing  vortices,  this  is  the  dominant  mechanism  for 
decrease  in  velocities  of  swirl  in  the  core  and  corresponding  growth  of  the 
core.  Axial  velocity  difference  between  the  core  of  the  trailing-vortex  and  i 
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the  surrounding  region  is  necessary  for  the  sustenance  of  turbulence  in  the 
^ vortex  core.  A theory  of  turbulent  trai ling-vortex  is  developed  on  the  basis 
of  these  mechanisms  and  the  results  are  compared  with  experimental  observation 
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